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ABSTRACT 

The redshifted 21cm signal of neutral hydrogen from the epoch of reionization (EoR) is extremely 
weak and its first detection is therefore expected to be statistical with first-generation low-frequency 
radio interferometers. In this letter we propose a method to extract the angular power spectrum of 
EoR from the visibility correlation coefficients pij(u,v), instead of the visibilities Vij(u, v) measured 
directly by radio interferometers in conventional algorithm. The visibility correlation coefficients are 
defined as Pij(u, v) = Vij(u, v)/ y/\Vu\\Vjj | by introducing the auto-correlation terms Vu and Vjj such 
that the angular power spectrum Ci can be obtained through Cg = 47t 2 Tq (\pij(u, v)\ 2 }, independently 
of the primary beams of antennas. This also removes partially the influence of receiver gains in the 
measurement of Cg because the amplitudes of the gains cancel each other out in the statistical average 
operation of (\pij{u, v)\ 2 }. The average brightness temperature To of extragalactic sources is used as 
the calibrator of Cg. Finally we demonstrate the feasibility of the novel method using the simulated 
sky maps as targets and the 21 CentiMeter Array (21CMA) as interferometer. 

Subject headings: cosmology: theory — diffuse radiation — intergalactic medium — methods: data 
analysis — techniques: interferometric 



1. INTRODUCTION 

While the redshifted 21cm emission/ absorption of neu- 
tral hydrogen provides a unique cosmological probe of the 
epoch of reionization (EoR) -, the last frontier of obser- 
vational cosmology, there are two primary challenges of 
measuring the EoR signatures with most of the dedicated 
radio facilities (e.g., 21CMA, LOFAR, LWA, MWA, PA- 
PER, SKA, etc.): First, the cosmic signal from EoR 
is deeply buried under the extremely bright foreground 
dominated by our Galaxy, extragalactic sources and tele- 
scope noise, and an unprecedented level of foreground 
removals down to five orders of magnitude should be re- 
quired in order to detect the cosmic signal (e.g., Madau 
et al. 1997; Zaldarriaga et al. 2004). Second, low- 
frequency radio interferometric measurements in current 
21cm experiments suffer from various instrumental con- 
taminations in addition to man-made radio- frequency in- 
terference. Outstanding among these are the frequency- 
dependent point spread function and field-of-view (also 
known as 'mode-mixing'), complexity of calibration, and 
bright source subtraction (e.g., Morales et al. 2006; Liu 
et al. 2009; Bowman et al. 2009; Datta et al. 2009,2010; 
Bernardi et al. 2010; Datta et al. 2010; Petrovic & 
Oh 2011). Yet, it is generally agreed among the 21cm 
cosmology community that the advent of many sophis- 
ticated techniques and algorithms in recent years helps 
overcoming these observational and technical hurdles, al- 
lowing to reach the desired detection sensitivity with the 
first generation of radio interferometers (for recent re- 
views see Furlanetto et al. 2006; Morales & Wyithe 2010; 
Pritchard & Loeb 2012; Zaroubi 2012). 



The theoretically predicted brightness temperature of 
21cm signal from EoR is only ~ 10 mK, and first cap- 
ture of such extremely weak signal is therefore expected 
to be statistical (Zaldarriaga et al. 2004). Most of the 
current 21cm experiments aiming to detect the EoR sig- 
nal are based on radio interferometric technique, which 
provides a direct measure of the Fourier component of 
the sky brightness I(s) convolved with the primary beam 
of the antennas -B«(s) towards direction s, often known 
as the visibility Vij(u) at a given baseline u in units of 
wavelength: 



Vy(u)=5* ff J / By(s)I(s)e 



— 27TIU-S j2 c 



(1) 



where gi and gj are the complex gain factors of antenna 
pair i and j , respectively. It can be easily shown that the 
angular power spectrum (Cg) of the sky brightness distri- 
bution I(s) can be constructed by the average value of the 
square of Vij (u) with the Fourier wavenumber £ = 2nu 
(White et al. 1999). In particular, under the assump- 
tion that the angular power spectrum varies rather slowly 
with scale relative to the Fourier component of the pri- 
mary beam Bij(u) = J-[Bij(s)], we obtain the commonly 
used formula in the estimation of the angular power spec- 
trum of low-frequency sky at a given frequency (Bharad- 
waj & Sethi 2001; Zaldarriaga et al. 2004; Bharadwaj & 
Ali 2005; Santos et al. 2005; Ali et al. 2008; Pen et al. 
2009; Paciga et al. 2011; Ghosh et al. 2011a,2011b): 
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However, most of the theoretical studies in literature im- 
plicitly assume that perfect calibrations are being made 
for radio interferometers, and gains and primary beam 
introduce no spectral and spatial structures to destroy 
the reconstruction of the EoR angular power spectrum. 



2 



Zheng et al. 



But in reality, to perform precise calibration in the ex- 
istence of ionospheric turbulence and to maintain high 
stability of radio instruments for rather a long integra- 
tion time arc extremely difficult for current 21cm exper- 
iments. This motivates us in this Letter to explore a 
possible remedy to overcome some of these shortcomings 
by using the visibility correlation coefficients instead of 
the visibilities. With the novel method it will be possible 
to statistically extract the EoR angular power spectrum 
from radio interferometric measurements, independently 
of the primary beam of antennas. The method also al- 
lows us to remove partially the influence of receiver gains 
in the statistical measurement of Ci because the ampli- 
tudes of the gains cancel each other out. 

2. FORMALISM 

We begin with the auto-correlation of voltages mea- 
sured at each antenna, which corresponds to the visibility 
with baseline u = 0: 



\Vu\ 2 = \g*\ 4 lt J d 2 s\B u (s)\ 



(3) 



This formula applies for discrete extragalactic sources 
only, and Jo denotes the mean brightness of the point 
sources rather than the total power. The latter is domi- 
nated by the Milky Way in our interested frequency range 
below 200 MHz. Now we define the visibility correlation 
coefficient such that 



(4) 



With this definition the amplitudes of the complex gains 
cancel each other out but the signal coherence remains. 
Following eq.(2) we take the average value of the square 

of Pij(u) 
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Using the 'frequency shift' theorem of the Fourier 
transform, we can rewrite the Fourier component of 
the primary beam Bij(u — u') as Bij(u — u') = 
J r [Bij(-s)e~ 2 ' Kiu ' s ]. According to the Parseval's theo- 
rem, the total power in the s domain or the u' domain 
should be the same. This yields 

/ dV|B«(u-u')| 2 = (6) 

J d 2 s\B l3 (-s)e- 2mu - s \ 2 = J d 2 s\B l3 (-s)\ 2 . 

The most crucial point of this equation, however, is that 
the power J d 2 u'\Bij(u — u')| 2 is actually independent 
of u. Namely, in the conventional estimate of angu- 
lar power spectrum from eq.(2) the primary beam term 
J d 2 u'\Bij(u — u')| 2 does not alter the shape of the 
power spectrum. This arises, of course, from the pre- 
sumption that the primary beam varies rather slowly 
with angular scale relative to the cosmic signal, a nec- 
essary prerequisite for taking out from the integral the 
Ct = 2-Ku term in eq.(2). Furthermore, the symmetric 
feature of the primary beam for interferometric array 



element in all the current 21cm experiments suggests 
that Bij(s) can be actually treated as an even func- 
tion, implying Bij(— s) = -By (s). Now, we apply the 
Cauthy-Schwarz inequality to the integral J d 2 s\Bij(s)\ 2 
by noticing By(s) = B;(s)5*(s), |B,(s)| 2 = |B«(s)| and 
|Bj(s)| 2 = \Bjj(s)\, and obtain 
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The equals sign holds when all antenna elements have 
identical primary beam J>i(s) = Bj(s). Finally, replacing 

the integral J ePu'|J?y(u — u')| 2 in eq.(5) by eqs.(6) and 
(7) and using again the Parseval's theorem we find 
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This allows us to evaluate the angular power spectrum 
C(, relative to the mean brightness J of point sources 
simply through 

C /=2wu =47T 2 jg<|py(u)| 2 ). (9) 

Consequently, we may use the average sky brightness of 
extragalactic sources over the primary beam of the an- 
tenna as the reference or calibration. While many so- 
phisticated models have been constructed for low fre- 
quency sky for the purpose of foreground removals in 
21cm experiments (e.g. de Oliveira-Costa et al. 2008; 
Jelic et al. 2008; Wilman et al. 2008; Vernstrom et 
al. 2011; Kogut 2012), for our purpose a simple power- 
law as the global sky brightness temperature model 
of extragalactic sources should suffice to calibrate the 
angular power spectrum at each frequency v. Tq — 

13.7K ( i50MHz ) • Thi s * s * ne best-fit mean sky tem- 
perature over a field of view of 20° x 20° at the frequency 
range 100-200 MHz based on the numerical simulations 
by Wilman et al (2008) after bright sources with fluxes 
greater than 10 mJy at 150 MHz are removed (see Section 
3). Note that Ci will represent the angular power spec- 
trum of the sky brightness temperature distribution if Jo 
is replaced by T in eq.(9): Ce=2nu = 47r 2 T 2 (|py(u)| 2 ). 

3. SIMULATION TEST 

In order to test the feasibility of extracting angular 
power spectrum of low-frequency sky with the novel 
method, we use the simulated sky maps by Wilman 
et al. (2008) as the input and the 21 CentiMeter Ar- 
ray (21CMA) as the radio interferometer. Following the 
same algorithm of Zheng et al. (2012), we generate a set 
of sky maps of 20° x 20° for different frequencies, con- 
taining five distinct radio source types out to redshift of 
z=20. We estimate the flux of each radio source at differ- 
ent observing frequencies between 100 MHz and 200 MHz 
by a running power-law in frequency: S cx i/ a + Aa 1 &( V / V *) ^ 
in which z/* is the characteristic frequency and is taken 
to be 150 MHz in this study. Finally, we exclude all the 
bright sources with fluxes Sisomhz exceeding 10 mJy at 
150 MHz which are assumed to be resolvable by current 
radio interferometers and can therefore be excised with 
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existing algorithms. We fix a 4096 2 grid for each of the 
simulated image, which gives rise to an angular resolu- 
tion of .3 (see Figure 1). 

21CMA, sited in Western China, is a ground-based 
meter-wave array designed to probe the EoR operating 
at frequencies from 50 MHz to 200 MHz. The array 
consists of 80 pods (or stations) with 127 log-period an- 
tennas for each, which are deployed in two perpendicular 
arms along east-west and north-south directions, respec- 
tively. Spacings between pods are integral multiples of 

20 m, with a maximum baseline of 2740 m along each 
baseline. The 21 CM A redundancy is being used for the 
purposes of not only calibrations but also statistical mea- 
surement of the angular power spectrum of EoR at spe- 
cific modes. In this work we choose the 40 pods of the 
east-west arm to generate the uv sampling towards the 
north celestial pole region, which reduces significantly 
the computing complexity in visibilities Vij{u,v) since 
only two-dimensional Fourier transform is involved for 
such a configuration. The corresponding baseline (i.e. uv 
sampling) distribution of the 780 pod pairs as interferom- 
eters for the 21CMA east-west arm is shown in Figure 2, 
among which there are only 127 independent baselines. 

The simulated sky map is convolved with the 21 CM A 
primary beam B,(s) which is identical for all the 
pods and takes approximately a Gaussian function with 
FWHM= 4°.26(i//100MHz)- 1 . We perform the Fourier 
transform of the simulated sky map modulated by (s) 
to produce the uv map in terms of eq.(l) by setting 
,9i = 9j = 1- The uv map is further sampled by the 

21 CM A east- west baselines shown in Figure 2. We now 
calculate the visibility correlation coefficients Pij(u, v) at 
each frequency channel for all the 780 baselines. This 
yields a set of Pij(u,v) measures at 127 independent 

Fourier modes I — 2tt\/u 2 + v 2 . Finally, the angular 
power spectra at these specific modes can be obtained 
using eq.(9), Ct=2nu = 47r 2 T ( 2 (\p tJ -,(u, v)| 2 ), in which the 
average sky brightness temperature for our case is de- 
termined by the unresolved extragalactic sources To = 
13.7K(f/150MHz) -2 - 74 because the Galactic foreground 
has not been included in the simulation and the bright 
sources have been already removed. 

Figure 3 shows the angular power spectra of the simu- 
lated sky maps at four frequencies ranging from 120 MHz 
to 180 MHz, together with the recovered ones at specific 
modes sampled by the 21CMA east-west baselines. The 
former are constructed directly using the Fourier trans- 
form of the simulated sky images without inclusion of 
any instrumental effects, while the latter are the recon- 
structed results from the 21CMA 'observations' based on 
the novel algorithm of eq.(9). It appears that the two re- 
sults show a remarkably good agreement. We have also 
demonstrated the measurement errors assuming an in- 
tegration time of 300 days, an observing bandwidth of 
1 MHz and the system noise of 300 K. In particular, 
the noise level at Ct has been suppressed by a factor of 
1/y/Ne for the redundant baseline of Nt equally spaced 
pods. Large error bars at small- and large-^ ends can be 
attributed to the arcminute-scale angular resolution due 
to the short baselines and the cosmic variance due to the 
small ficld-of-view of the 21 CM A, respectively. For com- 
parison the reconstructed angular power spectra directly 
from Figure 1 are also plotted in Figure 3. It appears 



that the Gaussian beam alters only the amplitude rather 
than the shape of the angular power spectrum in terms 
of eqs.(2) and (6). As is shown in Figure 3, such an 
amplitude effect has been corrected for when the visi- 
bility correlation coefficients pij(u,v) are used. Yet, in 
practice the primary beam can hardly be modeled by a 
perfect Gaussian function or other form of simple analyt- 
ical function. A careful calibration of the primary beam 
of antennas to a high degree of precision must be made. 
Imperfect and inaccurate calibrations of both spatial and 
spectral properties of the primary beam may lead to sig- 
nificant errors in reconstruction of the power spectrum of 
EoR for 21cm experiments. Employment of the visibility 
correlation coefficients in the statistical study of EoR al- 
lows us to eliminate concern about the calibration of the 
primary beam. While with our new algorithm we have 
successfully recovered the angular power spectra of the 
simulated radio foregrounds, the foregrounds should be 
eventually suppressed to the level below 10 mK, a min- 
imum requirement for extracting statistically the signa- 
tures of EoR. This can be achieved, for example, using 
the foreground removal technique suggested recently by 
Cho et al.(2012), which works straightforwardly with the 
angular power spectrum. We have tested the technique 
and found that the foregrounds can indeed be subtracted 
to the level below 10 mK. 

4. DISCUSSION AND CONCLUSIONS 

Instead of directly employing the visibilities in con- 
ventional interferometric measurements of the statistical 
fluctuations of EoR suggested in literature, we propose to 
work with the visibility correlation coefficients defined by 
Pij(u,v) = Vij(u,v)/y l /\Va\\Vjj\. This allows us to elim- 
inate the effect of primary beams of antennas and also 
partially reduce the influence of receiver gains on the sta- 
tistical extraction of the angular power spectrum of the 
low-frequency sky: Ce=2nu = 47T 2 T 2 (|p ii (w, v)\ 2 ). Yet, 
we need to calibrate the power spectrum using the av- 
erage brightness temperature T of extragalactic sources 
for low- frequencies below 200 MHz. Observationally, T 
has been determined so far to a degree of satisfaction at 
least for our purpose (e.g., Di Mattco ct al. 2002). 

Introduction of the visibility correlation coefficient 
Pij(u,v) does not change the coherence of original sig- 
nal. Hence the phase correction such as self-calibration 
and/or redundant calibration should be still made before 
combining Pij(u, v) data. Furthermore, bright sources in 
the field-of-view have to be removed to reduce the Pois- 
son noise in computation of the angular power spectrum. 
This also implies that the sidelobes of bright sources 
still remain as troublesome for reconstruction of angular 
power spectrum no matter whether Vij(u,v) or pij(u,v) 
is used. Another reason that bright sources should be 
excised before recovery of the angular power spectrum 
is the requirement of uniformality assumption for taking 
out from the integral the angular power spectrum term 
Ct in cq.(2). 

We have tested the feasibility of the novel method us- 
ing the simulated sky maps of 20° x 20° for extragalactic 
sources in low-frequencies as the targets and 21CMA as 
the radio interferometer. We have successfully recovered 
the angular power spectra of the foregrounds at specific 
Fourier modes sampled by the 127 independent baselines 
of the 21CMA east-west arm, after the bright sources 
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with fluxes of Sisomhz > 10 mJy are removed. While we 
have not included the gain fluctuations in sampling of 
the visibilities, the new method does allow us to remove 
the effect of the spatial response of the 21CMA anten- 
nas on the reconstruction of the angular power spectrum 
of the low-frequency sky. In combination of various so- 
phisticated foreground removal techniques developed in 
recent years especially in the power spectrum domain 
(e.g. Cho et al. 2012), we should be able to subtract 
the foreground to the level for statistical detection of the 



21cm signal from EoR. 
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Fig. 1. — Simulated sky map of 20° X 20° at v = 100 MHz convolved with the 21CMA antenna primary beam. All the bright sources 
SisomHz > 10 mjy have been removed. 
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Fig. 2. — Baseline distribution of the 780 pod pairs as interferometers in the 21CMA east-west arm, which also illustrates the uv sampling 
density when the horizontal axis is scaled in units of observing wavelength. Note that there are only 127 independent baselines among the 
780 pairs. 
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Fig. 3. — Angular power spectra of the simulated sky maps at four frequencies ranging from 120 MHz to 180 MHz, represented by 
ST = [l(2t + l)Q/47r] 1 /2. Solid lines are the results derived from the simulated maps without any observational and instrumental effects; 
blue circles are the reconstructed ST, in terms of the visibility correlation coefficients, measured at the 21CMA east-west independent 
baselines. For comparison the angular power spectra constructed directly from Fig.l are also shown (red circles), in which the Gaussian 
beam results in a decrease of the amplitude but does not alter the overall shape of ST. 



